Proposal of a one-dimensional electron gas in the steps at the LaA103-SrTi03 interface 
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The two-dimensional electron gas (2DEG) at the interface between LaAlOs (LAO) and SrTiOa 
(STO) has become one of the most fascinating and highly-debated oxide systems of recent times. 
Here we propose that a one-dimensional electron gas (IDEG) can be engineered at the step edges 
of the LAO/STO interface. These predictions are supported by first principles calculations and 
electrostatic modeling which elucidate the origin of the IDEG as an electronic reconstruction to 
compensate a net surface charge in the step edge. The results suggest a novel route to increasing 
the functional density in these electronic interfaces. 



The two-dimensional electron gas (2DEG) at the in- 
terface between two band insulators, LaAlOa (LAO) 
and SrTiOs (STO), has become the prototypical system 
for the study of novel electronic phases at oxide inter- 
faces. Since the original discovery [T] it has been shown 
that at low temperatures the 2DEG may become mag- 
netic [2], superconducting [3 , or both coexisting in one 
phase [HE]. At room temperature the LAO-STO inter- 
face may find applications in field-effect devices [6l [7], 
sensors [8 or nano photo-detectors [9^. 

The origin of the 2DEG is still under debate, but one 
potential mechanism, called the "electronic reconstruc- 
tion" [To], is currently the most popular [11]. At the 
heart of this mechanism is the notion that LAO has a 
non-zero polarization, P [12] [13], which for thin films, 
translates to having polar surfaces/interfaces along cer- 
tain crystallographic directions [TT, (including (001) - 
the most common growth direction in practice). Polar 
surfaces are electrostatically unstable but one possible 
route to charge compensate is via a transfer of electrons 
between the surface of the film and the interface. For a 
LAO film grown on a Ti02 terminated STO (001) sub- 
strate, complete compensation would amount to precisely 
0.5 e/ S (where S is the area of the (001) plane associated 
to one formula unit) transferring from the LAO surface to 
the interface, creating the 2DEG. Incomplete, or no com- 
pensation may occur for ultrathin LAO films, where the 
internal potential drop is lower than the relevant effective 
gap (in this case the band gap of STO). This is consistent 
with the observation of a critical LAO thickness, dc^ for 
2DEG appearance [6j . An alternative charge compensat- 
ing mechanism is via surface oxygen vacancies ^ which 
is also consistent with the observed dc (see Ref. fTS]). 

In this paper we propose that steps at the interface, 
defects which are observed in practically all epitaxial 
films (due to a small miscut of the substrate by angle 



^), can be utilized to produce one-dimensional electron 
gases (IDEG) residing at the step edge. We use a sim- 
ple model to argue that a step of one unit cell in height 
alters the polarity of the interface and surface of LAO, 
since the step edge itself has a net charge of e/2S (i.e. 
a charge line density of se/2S^ with s the step height as 
in Fig. 1). An equivalent, and more rigorous, picture 
takes the component perpendicular to the stepped sur- 
face of the formal polarization of LAO (Fig. 1), which 
is precisely (1,1,1) e/2S as shown later in the text. The 
electrostatics of the interface are studied with a simple 
capacitor plate model, supported with density functional 
theory (DFT) calculations, that shows an alteration of 
the interface carrier density with miscut angle. A IDEG 
is predicted for small miscut angles when the LAO thick- 
ness is just below dc{0 = 0), the critical thickness for 
metallicity in the zero miscut system. 

We begin with a formal consideration of the polarity 




FIG. 1: Schematic diagram illustrating the effect of steps on 
the net charge of the interface. Top: Polarization argument 
considering LAO has a polarization, Pq along the [100], [010] 
and [001] directions, with P the net of the stepped interface. 
Bottom: an equivalent picture with charge densities ±ac on 
the terrace and step edge. 
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FIG. 2: (Color online) Left: The DFT simulation cell of the 
LAO/STO superlattice. The two interface terrace termina- 
tions are both Ti02-LaO (n-type) whereas the interface step 
termination is n-type along the bottom interface and p-type 
(SrO-A102) along the top interface, (grey: oxygen, dark blue: 
Al, light blue: La, yellow: Ti, green: Sr) Center: (001) Lay- 
ered density of states near the Fermi level of the stepped su- 
perlattice (solid lines) and pristine superlattice (dashed line). 
Right: The change in the layer density of free electrons on 
introduction of steps in the superlattice. 



of stepped LAO surfaces. We define the miscut angle 
as tan6> = s/t, where s and t are the step and terrace 
distances respectively (see Fig. 1). The steps can be in- 
tentionally miscut along a certain direction [TBI [17]. We 
consider the [100] case here, but the arguments can be 
trivially applied to any miscut orientation. Atomic force 
microscopy (AFM) measurements of epitaxial LAO-STO 
suggest steps of 1 unit cell in height [161 Hi] • Neighboring 
terraces of the surface are then of the same termination 
(BO2 here), but the step edge may be terminated either 
AO or BO2 (non-polar vicinal LAO surfaces have been 
considered in Ref.lS). Perhaps the conceptually simplest 
picture is to take what has been called "chemical" [19] 
(or "compositional" [20]) charges, dc = O.Se/^, along the 
terrace and step edges (Fig. 1 bottom). The total net 
charge of the surface, ^c, is then ac{cos6±sm6). The ± 
depends on the step termination: addition if the step 
and terrace have the same terminations, and subtraction 
if they have opposite terminations. This result can be 
confirmed within the framework of the modern theory of 
polarization [2]- Taking the formal polarization, P, de- 
fined by [22 the position of the ion cores and mapping 
the electronic wavefunction on to Wannier centers [23 , 
one finds for cubic LAO P = V3 e/2S along the [111] di- 
rection [13. The component of this polarization normal 
to the stepped surface (Fig. 1 top), P, is the net sur- 



face charge density and precisely the same as the above 
"chemical" charge. For tetragonal LAO (as in the case of 
epitaxial thin films), a trivial correction may be applied 
to the above through a redefinition of S. 

To support this picture, density functional theory 
(DFT) calculations were performed on two model LAO- 
STO systems: n-n superlattices, and n-p films, n and 
p denote the termination of LAO terraces (n - LaO, p 
- AIO2). Note that in both cases, the steps are created 
by simply altering one of the supercell vectors - there 
is no change in stoichiometry through the introduction 
of steps. The DFT calculations were performed using 
the Wu-Cohen exchange correlation functional [24 as im- 
plemented in the Siesta code [25), Details of the 
pseudopotentials and atomic-orbitals basis are given in 
Ref. IT9I We begin with n-n superlattices (Fig. 2). 

Fig. 2 shows the effect of steps on the electronic struc- 
ture of the n-n superlattice. The superlattice was con- 
structed from 4.5 unit cells of LAO and 8.5 unit cells 
of STO stacked along the [001] direction. Steps were 
introduced by taking a supercell and altering the a su- 
percell vector to ta+sc (see Fig. 2 left). This produces 
n-type steps along one interface and p-type steps along 
the other while maintaining exactly the same stoichiom- 
etry. Fig. 2 shows the case of t=3 and s=l. In the 
system without steps (which we call the pristine system) 
the Ti 3d conduction electrons distribute evenly between 
the two n-type terraces - 0.5 e/S at each interface with 
a decay through the STO slab (Fig. 2 center), consistent 
with previous DFT calculations [271 HH] • The inclusion 
of steps is found to transfer electrons from one inter- 
face to the other (Fig. 2 right), accumulating near the 
n-type steps, and depleting near the p-type steps, with 
respect to the pristine system. The total density of elec- 
trons transferred is approximately O.le/S (20% of dc). 
This is consistent with the number required to screen the 
bound charge density of each interface, including the step 
charges, as predicted by the model. This model system 
clearly shows a large effect of steps on the interface free 
electron density and distribution. Note that there has 
been no change in stoichiometry nor is there any other 
mechanism suggesting an additional charge transfer. 

We now consider a system of more relevance for most 
experimentally grown systems in reality - n-p LAO films 
on STO substrates. The system consists of a slab ge- 
ometry of d LAO unit cells stacked on 5 unit cells of 
STO along the [001] direction including a 15 A vacuum 
layer. A dipole correction was included in the vacuum 
layer to simulate open-circuit boundary conditions, pre- 
venting unphysical macroscopic electric fields due to the 
asymmetry of the slab and the periodic boundary con- 
ditions. The introduction of steps was the same as in 
the above superlattice case, and chosen such that n-type 
steps reside in the n-type terrace, and p-type steps reside 
in the p-type terrace. We begin with the case of d = 3 
unit cells, which was found to be beyond the pristine 
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FIG. 3: (Color online) Top: Interface free electron density, 
cr, as a function of the miscut angle, ^, for a LAO (3 unit 
cells) /STO film geometry. Points represent densities calcu- 
lated from DFT, and lines are electrostatic modeling under 
various assumptions (see text). Middle: Predicted schematic 
electronic phase diagram at the n-type (step and terrace) ter- 
minated LAO/STO interface as a function of the miscut angle, 
^, and the LAO thickness, d. Bottom: An isosurface of the 
Ti electron density (yellow clouds) at an n-type stepped 
LAO/STO interface. Blue and pink cages represent oxygen 
octahedra around Ti and Al atoms respectively. 



critical thickness within DFT for the onset of metallicity, 
and then alter the miscut angle, 6>, through varying t (5 
is again fixed at 1 unit cell to resemble that seen in ex- 
periments). The total interface free electron density, a, 
is calculated from the projected density of states, which 
allows one to easily separate free and bound charge due 
to the existence of a gap. The results are depicted as 



points in Fig. 3 top, which shows a clear variation of 
free charge density with miscut angle, going through a 
maximum (nearly 200% of = 0)) at intermediate an- 
gles. Three curves are compared in Fig. 3 top, each with 
various model approximations, which we discuss next. 

To model <j as a function of d we use (i) Gauss' Law, 
(Jc — cr = e£, where e is the LAO dielectric constant and 
E the LAO internal electric field, and {%%) voltage drop 
pinning, 5= A/d [19 (A is the STO band gap), obtaining 



eA 

d 



(1) 



with Gc = crc(cos^+sin^) the net charge of the 
stepped interface as discussed earlier, and d = d{l -\- 
Ctan^)/(l+tan6>), the average shortest distance between 
the stepped plates {C is a geometrical constant depend- 
ing only on the ratio s/d). Eq. 1 produces the bold line 
curve in Fig. 3 top, which is in good agreement with 
the DFT calculations. Importantly, both the effect of 
the steps on the net polarization, ac, and on the effective 
LAO thickness, d, is required to agree with the DFT cal- 
culations. Substituting simply the terrace net charge dc 
and or the LAO thickness d results in disagreement with 
DFT (dotted and thin hues Fig. 3 top). 

With strong support from the DFT calculations, we 
explore the model further. Firstly the critical thickness 
for metallicity as a function of miscut angle, dc{0). This 
is determined from Eq. 1 by setting a to zero and rear- 
ranging for d. A schematic of this transition is shown in 
Fig. 3 middle, for small ^, as the curved line separating 
the insulating phase. The degree of curvature away from 
horizontal has been exaggerated for clarity. In fact, using 
the experimental values for e and A and cic = O.Se/^, it 
reaches a minimum of about 97% of dc{0 = 0) ai ^ 17° . 
Interestingly, this minimum can be altered through epi- 
taxial strain since the step edge charge density relies on 
the out-of-plane lattice parameter. Fig. 3 middle shows 
three other phases above this transition curve - 2DEG, 
anisotropic 2DEG and IDEG. Above the pristine crit- 
ical thickness (horizontal line) an anisotropic 2DEG is 
predicted at any non-zero miscut angle. This is due to 
the non-isotropic net charge density of the interface (see 
later for support from DFT). Between the pristine critical 
thickness and dc{0)^ and for small ^, a IDEG phase is pre- 
dicted. Introducing steps just below the pristine critical 
thickness produces carriers in the step edge, to screen the 
excess interface bound charge. In the limit of very small 
these are effectively isolated steps. Considering one sin- 
gle step, the excess chemical charge is located around the 
step for the same reason that the chemical charge associ- 
ated to the pristine LAO/STO (001) interface is located 
close to the interface. Since the free carrier charge being 
transferred is less than the excess chemical charge, the 
potential well generated by the latter confines the carri- 
ers within some width around the step in the direction 
normal to the step, thereby defining a IDEG. A phase 
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transition with increasing miscut angle occurs when the 
IDEGs, with a finite width, begin to overlap with next 
neighbor steps. At this point, the IDEGs are electroni- 
cally indistinguishable from the anisotropic 2DEG. 

In an attempt to observe the IDEG, DFT calculations 
were performed on the n—p film system (see above), with 

as 2 unit cells, which is lower than the pristine critical 
thickness, and t=5, s=l. Fig. 3 bottom shows an isosur- 
face of the Ti 3d electron density. The electron density 
is clearly inhomogeneous, and accumulates near the step 
edge. The density decays from the step edge in to the 
terrace, but not to zero. The phase is the anisotropic 
2DEG of Fig. 3 middle, below the dotted line. By re- 
ducing the miscut angle in the simulation, it should be 
possible to find the IDEG, however we are limited to not 
too small miscut angles due to system size. 

Experimental studies of stepped LAO-STO interfaces 
are few, but we highlight two here. Anisotropic elec- 
trical transport has been observed, but explained by 
anisotropic carrier mobility, not density [16^. Interest- 
ingly, the strongest anisotropy was observed at lowest 
miscut angles, in disagreement with their step scattering 
model. Elsewhere, steps were found to affect the carrier 
mobility and density of the 2DEG The carrier den- 
sity was found to decrease with miscut angle, suggesting 
the possibility of p-type step formation. 

We hope this work will stimulate the study of stepped 
surfaces and the search for IDEGs. Likely candidate sys- 
tems include; {i) n-type step and terrace terminations, 
{ii) a theoretical pristine critical thickness of only just 
beyond an integer number of LAO unit cells, {in) a max- 
imised curvature of dc{0) and {iv) a small miscut an- 
gle to disconnect the IDEGs from successive steps, but 
not too small since the effect of the step polarization is 
then minimized. In practice this may mean; {i) surface 
chemical treatment to alter the step edge termination, 
{ii) engineering the pristine critical thickness (through 
for example strain [19^, ^29^, electric field [6l [7j, mate- 
rial [30] or atmospheric environment if the electronic re- 
construction instead occurs via surface chemical redox 
reactions ^ llSj), and {Hi) increasing the step bound 
charge through epitaxial strain, and miscut orientation. 

In conclusion we have studied the effects of steps at 
the LAO-STO interface on the electronic structure. A 
simple model of polarization and electrostatics is in good 
agreement with DFT calculations. The model predicts 
the existence of a IDEG in the step edge at the interface. 
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